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Te covered Si(OOl): a variable surface reconstruction 
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At a given temperature, clean and adatom covered silicon surfaces usually exhibit well-defined reconstruction 
patterns. Our finite temperature ab-initio molecular dynamics calculations show that the tellurium covered 
Si(OOl) surface is an exception. Soft longitudinal modes of surface phonons due to the strongly anharmonic 
potential of the bridged tellurium atoms prevent the reconstruction structure from attaining any permanent, 
two dimensional periodic geometry. This explains why experiments attempting to find a definite model for 
the reconstruction have reached conflicting conclusions. 
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Relaxation and reconstruction of clean and adatom 
covered surfaces is an active field of study. Tremendous 
efforts have been devoted to observing and understand- 
ing how the symmetry and atomic configurations of sur- 
faces change, and how these changes affect the chemical 
and physical properties of surfaces. Particular atomic 
structures with well-defined reconstruction geometry are 
verified and sometimes predicted, by performing static 
total energy calculations at T=0 K. Atomic configura- 
tions corresponding to the global or local minima on a 
Born-Oppenheimer surface are then attributed to stable 
surface structures. 

In an effort to promote technology by growing crys- 
tals with the minimum possible defects, the clean and 
adatom covered surfaces of silicon have been thoroughly 
investigated. Recent studies on the growth of silicon sur- 
faces have shown that atoms like As, Sb, Te are good 
surfactants preventing island formation and hence aiding 
the layer by layer growth. |ffl-H In addition, the goal of 
combining large infra-red detector arrays with relatively 
cheap and well developed Si integrated circuit technology 
has led to the growth of HgCdTe on Si(001) surface. Be- 
cause of a large (~ 19%) lattice mismatch, thick buffers 
layers of CdTe must be grown before growing active lay- 
ers of HgCdTe. The atomic configuration of adsorbed 
Te, which forms the first layer grown on the bare Si(001) 
surface, is crucial for the fabrication of high performance 
devices, since it determines how growth nucleates and 
hence the quality of the epilayer. || 

The adsorption of Te on Si(001) surfaces for different 
coverages and the resulting atomic geometries have been 
studied by different surface techniques. While the ener- 
getically favorable adsorption sites are well understood 
at low Te coverage, |?],|9| the models of surface recon- 
structions proposed for near monolayer coverage (O ~1) 
have been at variance. [^o|-^2| Most of the LEED experi- 
ments observe a (lxl) structure up to temperatures high 
enough to desorb Te from the surface. jlO| Other exper- 
iments report different structures. For example, Tamiya 
et al. |1 1| found the transition from the low tempera- 



Wiame et al. [O observed a (2x1) symmetry in their 
STM images even at room temperature and proposed an 
atomistic model involving Te-Te dimers. The differing 
of the reconstruction geometries from one experiment to 
another is a puzzling and an uncommon situation. 

In this work, we explain this puzzling situation us- 
ing the results of finite temperature ab-initio molecular 
dynamics (MD) calculations. We investigate the ener- 
getics of Te adsorption starting from very low coverage 
(O =0.0625) up to a monolayer coverage (0 =1). We 
first determine the binding energies of a single Te atom 
adsorbed at the special sites of the unit cell for <C 1. 
We describe how the original Si-Si dimer bonds of the 
Si(001)-(2x 1) surface are broken, and how eventually the 
surface is covered by Te atoms. We also examine the 
possibility of two adjacent adsorbed Te atoms forming a 
dimer bond to give a (2x1) reconstruction. Other possi- 
ble higher order reconstruction geometries are searched 
by a finite temperature, ab initio molecular dynamics 
method. We find that uncorrelated lateral excursions of 
bridged Te atoms in flat potential wells hinder the obser- 
vation of any definitive surface reconstruction pattern at 
finite temperatures. 

Calculations were carried out within the density func- 
tional approach using the Vienna Ab-initio Simulation 
Package (VASP) (l|. The wave functions are expressed 
by plane waves with the cutoff energy |k + G| 2 < 250 eV. 
The Brillouin zone (BZ) integration is performed by using 
the Monkhorst-Pack scheme |]| with (2x2x1), (2x8x1) 
and (4x8x1) special points for (4x4), (4x1) and (2x1) 
cells, respectively. The convergence with respect to the 
energy cutoff and number of k-points was tested. Ionic 
potentials are represented by ultra-soft Vanderbilt type 
pseudopotentials |L6| and results are obtained within gen- 
eralized gradient approximation |l7| for a fully relaxed 
atomic structure. The preconditioned conjugate gradient 
method is used for wave function optimization and the 
conjugate gradient method for ionic relaxation at T =0 
K. At finite temperatures, the Nose-Hoover thermostat 
fl8f is employed for constant temperature dynamics of 



The time step in MD calculations, At is chosen such that 
typical phonon time period is divided into a few tens of 
time steps. We picked At to be 2 fs to ensure that the 
ionic trajectories are smooth. 

The Si(OOl) surface is represented by a repeating slab 
geometry. Each slab contains 5 Si(OOl) atomic planes and 
hydrogen atoms passivating the Si atoms at the bottom 
of the slab. Consecutive slabs are separated by a vacuum 
space of 9 A. For calculations at T =0 K, Si atoms in the 
top four atomic layers are allowed to relax, while the 
bottom Si atoms and passivating hydrogens are fixed to 
simulate bulk-like termination. 

In finite temperature calculations, all atoms, including 
Si and H atoms in the bottom layer, are allowed to move 
to avoid a large temperature gradient. Lattice param- 
eters are expanded according to the temperature under 
study using the experimental thermal expansion coeffi- 
cient in order to prevent the lattice from experiencing 
internal thermal strain. We reproduced the energetics 
and geometry of the c(4 x 2), p(2 x 2) and p(2 x 1) re- 
constructions of a clean Si(001) surface using the above 
parameters. fl4|] 

The binding energy of a single Te adsorbed on the spe- 
cial (on-top T, cave C, hollow H, and bridge B) sites on 
the clean Si(001) surface are calculated using a supercell 
consisting of eight (2x1) cells. The large size of the su- 
percell ensures that the interaction between the adsorbed 
Te atoms is negligible so that results can represent low Te 
coverage. In Fig. la, only one (2x1) cell of the supercell 
is shown. The binding energies are found to be T: 4.5 eV, 
C: 3.5 eV, H: 3.4 eV, and B: 3.2 eV. These binding en- 
ergies were calculated for fully relaxed structures at T=0 
K. Apparently, the most energetic site at low coverage is 
the on-top site, where a Te atom above the dimer bond of 
the clean Si(001)-(2xl) surface is bonded to two Si atoms 
of the same dimer bond. This is consistent with our in- 
tuitive chemical notion that Te(5p 4 ) tries to fill its outer- 
most p-shell by coordinating with two surface Si atoms. 
Our result is also in agreement with STM images. Q By 
considering only two special sites, Takeuchi |19|] found 
the on-top site to be energetically more favorable than 
the bridge sites by 0.8 eV. We examined the stability of 
the Te atom adsorbed at the on-top site for higher cov- 
erages. For O =0.5, Te atoms adsorbed 2.25 A above 
each surface dimer bond were found stable, except that 
the underlying Si-Si dimer bond is elongated marginally 
and the dimer asymmetry is removed. The Si-Te bond 
length is 2.53 A which is close to the sum of the Si and 
Te covalent radii and in excellent agreement with exper- 
iment. (20) 

A monolayer coverage of Te (i.e. O =1) is the most 
critical insofar as the controversy regarding the surface 
reconstruction is concerned. We attempt to resolve the 
controversy by addressing the following issues which are 
not settled yet. These are: i) How does the atomic con- 
figuration of the surface change with increasing >0.5 ? 
ii) Can two adjacent Te atoms on the surface dimerize at 



struction and how does the surface structure vary with 
temperature at ~ 1? To address the first question, we 
begin with an initial configuration where one Te is ad- 
sorbed at the T-site and the second one at the B site on 
the Si(001)-(2x 1) surface, and let this structure relax at 
T=0 K. The occupation of the B site at high Te cover- 
age is consistent with experiments. @j2lj In reaching the 
stable structure, the Te atoms form directional bonds 
with surface Si atoms while Si-Si dimer bonds elongate 
and eventually break. It appears that each Si-Si dimer 
bond is broken to form four new Si-Te bonds. In the fi- 
nal stable structure, Si atoms of the broken dimer bond 
are pushed to their bulk positions reforming the outer- 
most, bulklike Si(001) atomic plane. Each adsorbed Te 
atom is connected to the substrate with two Te-Si bonds 
of length 2.53 A. At the end, a metallic Te(001) atomic 
plane forms 1.65 A above the Si substrate with a bind- 
ing energy of 4.28 eV per Te atom relative to the clean 
Si(001)-(2x 1) surface and free Te atom. Figure lb de- 
scribes the atomic positions of this ideal (lxl) structure 
of the Te monolayer on the Si surface. The charge density 
contour plots in Fig. lc indicate that the bond is direc- 
tional. The maximum of the charge occurs between Si 
and Te, but closer to Te. |22| In spite of the directional 
Te-Si bonds, the surface of Te covered Si(001) surface 
is metallic with a finite density of states at the Fermi 
level. From force calculations we find that the Te atoms 
are robust against displacements along the [110] (or x-) 
direction in the plane of Si- Te-Si bonds. 

Our calculations for a free Te2 molecule predict a bind- 
ing energy of 4.41 eV and a bond length of 2.56 A. This 
suggests the possibility that two adjacent Te atoms on the 
Si(001) surface may experience energy benefit by form- 
ing a Te-Te dimer bond by moving towards each other in 
the y-direction. (See Fig. lb) Such a dimerization can 
once again lead to a (2x1) reconstruction. As a matter 
of fact, a similar As- As dimerization is known to occur 
on As covered Si(001) and Ge(001) surfaces. [|g| To 
test whether Te-Te dimerization can occur and to an- 
swer point (ii), an initial structure with a Te-Te distance 
of 3 A (which is greater than the bond length of Te2, but 
smaller than the undimerized distance in the (lxl) struc- 
ture) is relaxed at T=0 K. Upon relaxation, Te atoms 
moved away from each other so that the tilted Si- Te-Si 
plane became perpendicular to the surface and the to- 
tal energy of the system is lowered significantly. The 
analysis of the charge density in a (001) (or xy) plane 
passing through the Te atoms suggests that the forma- 
tion of strong Te-Si bonds excludes the bonding between 
two adjacent Te atoms (Fig. Id). Simple valence argu- 
ments also suggests that Te being divalent, would tend 
to avoid bonding with three other atoms. 

To address the most significant question (in) posed 
above, one must consider the reconstruction at =1 
which may involve complex and concerted rearrange- 
ments of the substrate and adsorbate atoms at high tem- 
perature. To access all possible reconstruction geometries 
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ysis at T=0 K, we performed a finite temperature, ab- 
initio molecular dynamics calculations at T=600 K and 
T=1000 K using (4x1) supercell geometry. Q Figure 
2 illustrates the displacements of Te atoms in a (4x1) 
supercell at T=600 K. The time variation of the mean 
squared planar displacements, (ttjj) = |Eti( u ii+ u j,i)> 
(w's are the displacements of the atoms from their ideal 
lattice positions) shows that the system is sufficiently 
thermalized within ~ 1 ps. We note that the displace- 
ments along the x-direction, u x ,i=iA(t), are small since 
the bridged Si-Te-Si bonds are robust. The average of 
the perpendicular positions of Te atoms on the surface, 
(z) = j YU=i z « W' an d a l so those of 8 hydrogen atoms 
at the bottom drifts along the z-direction with the same 
negative velocity, d(z)/dt ~-0.7 A/ps. In addition to 
this spurious translation of the unit cell, the displace- 
ment of each Te atom, u Zt i(t), oscillates with decreasing 
amplitude and without any correlation with the other Te 
atoms. 

The displacement along the [110] (or y-) direction, 
Uy t i(t), is large and can be relevant for a particular recon- 
struction structure. After the thermalization of the sys- 
tem, ity t i(t) becomes oscillatory and quasi periodic with 
periods of the order of ~ 1.0 ps. The behavior illustrated 
in Fig. 2 is reminiscent of the surface longitudinal acous- 
tic mode due to Te rows. The amplitudes of oscillations 
vary between 0.4 A to 0.7 A resulting in lateral excur- 
sions (as large as 1.4 A) of Te rows along [110] direction. 
To enhance the statistics, we performed the same calcu- 
lation at T=1000 K. The adsorbed Te atoms execute sim- 
ilar motions, only with larger amplitudes, at this higher 
temperature. 

These excursions or displacements of adjacent rows do 
not display any correlation. Moreover, they are time de- 
pendent. The random and uncorrelated nature of the 
displacements prevents us from deducing a well-defined 
reconstruction pattern. Such excursions of Te rows along 
the [110] direction would not give rise to any resolvable 
pattern in the LEED and STM images. For example, 
since the period of oscillations are much shorter than the 
characteristic scan time of STM, the STM images taken 
at finite temperature would indicate disordered (lxl) re- 
construction. 

For adsorbed Te rows to execute large amplitude excur- 
sions with low frequency at T=600 K is unusual and sug- 
gests rather soft and non-Hookian (nonlinear) force con- 
stants in this direction. In fact, as seen in Fig. 3, the total 
energy remains practically unchanged for a displacement 
of the Te rows of u y ~ ± 0.5 A. For the displacement of 
adjacent rows in opposite directions, Et(u v ) resembles a 
double well potential with a broad maximum at u y =0 
and a shallow minimum on cither sides. The barrier be- 
tween these two minima is very low, almost at the accu- 
racy limit of the present calculations (7 meV). This sug- 
gests that adjacent Te rows are displaced by ~0.25 A in 
opposite directions, forming a zigzag chain of Te atoms 
on the [110] direction and leading to a (2x1) surface re- 



construction at T=0 K. Interestingly, except for the dis- 
appearance of the weak double well form, the variation of 
the total energy with u v remains essentially unaltered if 
the adjacent Te rows are displaced in the same direction. 
This implies that, at finite temperatures, Te rows can eas- 
ily traverse the weak barrier and execute random (uncor- 
related) displacements. This situation is consistent with 
the results of finite temperature MD calculations sum- 
marized in Fig. 2. Since the potential energy well is so 
flat, the positions of Te atoms would be easily modified 
by the tip-sample interaction in STM experiments. The 
total energy curve in Fig. 3, is a fit to an analytical form 
E T (u y ) = cm* + [3u 4 y + ju y , (with a =0.3024 eV/A 2 , 
(3 =0.6242 eV/A 4 , 7 = -0.2087 eV/A 6 ) and reflects 
strong anharmonicity (nonlinearity in force constants) of 
the potential wells wherein Te atoms move. pj| 

In summary, we have found that Te atoms adsorb 
above the Si-Si dimer bonds at low coverage. There is 
no energy benefit for forming Te dimers at any coverage. 
At monolayer coverage, the potential wells for Te atoms 
are rather flat and strongly anharmonic along the [110] 
direction. There is almost no barrier for the Te rows on 
the surface to make significant excursions relative to their 
ideal positions along the [110] direction. First principle 
finite temperature calculations indicate that the displace- 
ments of Te rows are uncorrelated, lacking any definitive 
reconstruction pattern. 



[1] M. Copel, M. C. Reuter, E. Kaxiras, and R. M. Tromp, 

Phys. Rev. Lett. 63, 632 (1989). 
[2] M. Tromp and M. C. Reuter, Phys. Rev. Lett. 68, 954 

(1992). 

[3] S. Higuchi and Y. Nakanishi, Surf. Sci. 254, L465 (1991); 

H. J. Osten, J. Klatt, G. Lippert, E. Bugiel and S. 

Higuchi, J. Appl. Phys. 74, 2507 (1993). 
[4] N. Grandjean and J. Massies, Phys. Rev. B53, R13231 

(1996). 

[5] M. R. Bennett, C. J. Dunscombe, A. A. Cafollo, J. W. 
Cairns, J. E. Macdonald, and R. H. Williams, Surf. Sci. 
380, 178 (1997). 

[6] C. D. Consorte et. al. Phys. Rev. B63, R041301-1 (2001). 

[7] M. R. Bennett, A. A. Cafolla, J. W. Cairns, C. J. Dun- 
scombe, and R. H. Williams, Surf. Sci. 360, 187 (1996). 

[8] A. Almeida, Y. P. Chen, J. P. Faurie, S. Sivananthan, D. 
J. Smith and S. C. Tsen, J. Electron. Mater. 25, 1402 
(1996). 

[9] S. A. Yoshikawa, J. Nogami, C. F. Quate and P. Pianetta, 

Surf. Sci. 321, L183 (1994). 
[10] S. Di Nardo, L. Lozzi, M. Passacantando, P. Picozzi and 

S. Santucci, Surf. Sci. 331-333, 569 (1995). 
[11] K. Tamiya, T. Ohtani, Y. Takeda, T. Urano and S. 

Hongo, Surf. Sci. 408 (1998). 
[12] F. Wiame, G. Mahot, S. Sivananthan, S. Rujirawat, R. 



(1999). 

[13] G. Kresse and J. Hafner, Phys. Rev. B47, R558 (1993); 
G. Kresse and J. Furtmiiller, Phys. Rev. B54, 11169 
(1996). 

[14] A. Ramstad, G. Brocks, P. J. Kelly, Phys. Rev. B51, 
14504 (1995). 

[15] H. J. Monkhorst and J. D. Pack, Phys. Rev. B13, 5188 
(1976). 

[16] D. Vanderbilt, Phys. Rev. B41, 7892 (1990); G. Kresse 
and J. Hafner, J. Phys.: Condens. Matter, 6, 8245 (1994). 

[17] J. P. Perdew and Y. Wang, Phys. Rev. B46, 6671 (1992). 

[18] S. Nose, J. Chem. Phys. 81, 511 (1984). 

[19] N. Takeuchi, Phys. Rev. B60, 4796 (1999). 

[20] S. R. Burgess, B. C. C. Cowie, S. P. Wilks, P. R. Dunstan, 
C. J. Dunscombe and R. H. Williams, Appl. Surf. Sci. 
104-105, 152 (1996). 

[21] As compared to the C, and H site adsorptions, the B site 
adsorption is not energetically most favorable. However, 
the energetics change upon coadsorption with the T site. 

[22] According to Pauling's scale, Te is more electronegative 
than Si (i.e. \Ta =2.0 and \Si =1-8). Therefore, one ex- 
pects that charge is transferred from Si to Te. This is 
consistent with the atomic configurations and occupan- 
cies, Te(5s 2 5p 4 ) and Si(3s 2 3p 2 ). 

[23] R. I. G. Uhrberg, R. D. Bringas, R. Z. Bachrach and J. 
E. Northrup, Phys. Rev. Lett. 56, 520 (1986). 

[24] We recognize that high order reconstructions cannot be 
revealed by using a (4x1) supercell. For the Te covered 
Si(001) no high order reconstruction have been reported 
at e ~ 1. 

[25] Recently, it has been shown that strong anharmonicity 
of the potential energy is the cause of high temperature 
superconductivity in MgB2 crystals. T. Yildirim et al. 
Phys. Rev. Lett, (in press). 

FIG. 1. (a) The unit cell of the Si(001)-(2x 1) surface. 
The possible sites for the adsorption of Te at very low 
are marked by X. (b) The (lxl) structure of the Te cov- 
ered Si(001) surface, (c) Charge density contour plots of the 
Si-Te-Si bonds with arrows showing the direction of increas- 
ing charge density, (d) Charge density contour plots on the 
(001) plane passing through the Te atoms. Large filled, large 
empty, small empty, and smallest empty circles denote Te, 
first layer Si, second layer Si, and third layer Si atoms, re- 
spectively. The thick lines between circles indicate bonds. 
x-, {/-, and z-axis are parallel to the [110], [110], and [001] 
directions, respectively. The lattice constant a =3.84 A. 

FIG. 2. Time variation of the displacements of the Te 
atoms (u x , u y , and u z ) from their ideal lattice positions in a 
(4x1) supercell calculated at T=600K. The left panel shows 
the supercell. At t=0, all the atoms are at their ideal lattice 
positions. 



FIG. 3. Variation of total energy with the displacements of 
the Te row, u y , calculated at T=0 K. Each data point corre- 
sponds to a fully relaxed structure under a given displacement 
of the Te rows along the [110] direction. The thick line is for 
the adjacent rows moving in opposite directions forming a 
zigzag pattern. O's correspond to the Te rows moving in the 
same direction. An analytical fit, Et(u v ), to the total energy 
values depicted by the O's is shown by the broken line. Ener- 
gies are measured with respect to the perfect (lxl) surface 
for u v — 0. 
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